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,-_i ' Abstract 

0\ '. 

Qv^ , The paper deals with the one-dimensional parabolic potential barrier V{x) = 

(— I I Vo — ?n,7^x^/2, as a model of an unstable system in quantum mechanics. The time- 

Qh' independent Schrodinger equation for this model is set up as the eigenvalue problem 

(-H , in Gel'fand triplet and its exact solutions are expressed by generalized eigenfunc- 

"^ I tions belonging to complex energy eigenvalues Vq =f i^n whose imaginary parts are 

C^ ■ quantized as Fn = {n + l/2)fi'y. Under the assumption that time factors of an unsta- 

• • . ble system are square integrable, we provide a probabilistic interpretation of them. 

. !^ I This assumption leads to the separation of the domain of the time evolution, namely 

/\ • all the time factors belonging to the complex energy eigenvalues Vq — iFn exist on 

^ . the future part and all those belonging to the complex energy eigenvalues Vq + iFn 

exist on the past part. In this model the physical energy distributions worked out 

from these time factors are found to be the Breit-Wigner resonance formulas. The 

half-widths of these physical energy distributions are determined by the imaginary 

parts of complex energy eigenvalues, and hence they are also quantized. 



1 Introduction 

It is well known that time- dependent wave functions of non-relativistic quantum mechanics 
satisfy the time-dependent Schrodinger equation 

tn^^^P{t,q) = H^P{t,q), (1.1) 

where H is the Hamiltonian of the system expressed in the coordinate representation. Ac- 
cording to Born's probabilistic interpretation, the square of modulus of a time-dependent 
wave function ip^t, q) multiplied by a small volume Aq of the coordinate space, 

m,q)fAq (1.2) 

is the relative probability of the particle being within the small volume Aq of the point q 
at the time t. If the time- dependent wave function il){t, q) is normalized, 

m,q)fdq = l, (1.3) 

then the formula (|1 . 2| ) gives directly the probability. In this case the expectation value of 
any observable O is 

{0)t= U{t,qrdi^{t,q)dq, (1.4) 



which is a function of t. When a potential energy does not involve the time t explicitly, 
we can separate variables of a time-dependent wave function 

V'(t,g)=T(t)M(g). (1.5) 

For the state of a energy eigenvalue E, the time factor T(t) is of the form 

T{t) = Ae-'^*/^ (1.6) 

where A G C is independent of t, and equation ( |1.1| ) for it reduces to the time-independent 
Schrodinger equation 

Hu{q) = Eu{q). (1.7) 

In the case when the energy eigenvalue E of (|1.7|) is a real number, the formula (|1.2| ) will 
reduce to 

\u{q)fAq, (1.8) 



where u{q) contains A as a numerical coefficient. We then have, instead of ( |1.3|) 

\u{q)fdq = l, (1.9) 



and instead of ( |1.4| ) we get 

{0), = j u{qrdu{q)dq. (1.10) 

Thus the probabihty is independent of the time. We call a state belonging to the real 
energy eigenvalue a stationary state [|l], EL Q. If we cannot normalize an eigenfunction 
M(g), we obtain the expectation value of O, as the generalization of ( |1.10| ) 



(0),= / u{qrOu{q)dq / \u{q)\' dq. (^J^) 



However, the energy eigenvalue E can be a complex number for an unstable system 
in which the potential energy do not have a stable stationary point. The eigenfunction 
u{q) of the Hamiltonian H belonging to the complex energy eigenvalue is not expressible 
in terms of the vector of a Hilbert space, since H is a, self-adjoint operator on a Hilbert 
space. The problems of the mathematical formulations of continuous spectrum or complex 
eigenvalues has been treated by Bohm and Gadella |Q on the basis of a rigged Hilbert 
space or a Gel' f and triplet 0. They suggested by the analysis of S'-matrices that the 
states belonging to continuous spectrum or complex eigenvalues are expressible in terms 
of the generalized functions of a Gel'fand triplet, and their time evolution associates with 
a one-parameter semigroup. 

The object of the present paper is to show that the exact solution of the time- 



independent Schrodinger equation ( |1.7|) for an unstable potential energy becomes the 
generalized function belonging to the complex energy eigenvalue. As a model Hamilto- 
nian of the unstable system, we take the one-dimensional parabolic potential barrier. The 
collision problem of this model Hamiltonian was first solved by Kemble by the method 
of the WKB approximation. We consider the eigenvalue problem of this model Hamilto- 



nian in § 1^ at one instant of time. In § 2J^ this problem is worked out in the coordinate 



representation. We show that the energy eigenvalues are complex numbers and the corre- 
sponding eigenfunctions are expressible in terms of the generalized functions of a Gel'fand 
triplet. In § ^^ the same problem is dealt with in the momentum representation. 

Physically these complex eigenvalues play an important part in the resonance scatter- 
ing. Since a state of resonance is only approximately stationary, its mean lifetime will be 
finite. This means the probability of the particle being inside the system tends to zero 
after a sufficient lapse of time, and the generalized eigenfuction u{q) remains finite at 
|g| -^ oo. This is a physical reason why u{q) is not expressible in a Hilbert vector and 
does not satisfy the normalizing condition ( |1.9| ). Such a state is called a non- stationary 
state [0, ^ or a quasi- stationary state [^. Let us now try to assume the normalizing 
condition for the time factor 

\T{t)\^dt = l, (1.11) 



since the mean lifetime of the non- stationary state is finite. This condition will be seen 
later (see § p.2|) to be connected with the fact that the Breit-Wigner resonance formula is 
normalized. We can now proceed to introduce a physical interpretation of the time factor 
T(t) such that the absolute probability of the particle being during the time interval t to 
t + At is given by 

|T(t)|Mt. (1.12) 

This interpretation provides a physical description of a process of time for a non-stationary 
state, say from growth to decay. When this interpretation is adapted to Born's one, the 
square of modulus of ipit, q) corresponding to a non-stationary state, multiplied by AtAq, 

m,q)fAtAq (1.13) 

gives the relative probability of the particle being within the small volume Aq of the point 
q during the time interval t to t + At. The time average of the expectation value of O 
will be 

(O) = jjijit,qydi;{t,q)dtdq/j\uiq)fdq (1.14) 

as the generalization of ( |1.10f ). 

In § ^ the motion of the unstable system consisting of the one-dimensional parabolic 
potential barrier is considered. The above theory is applied to this motion in § p.l| . In 
§ p]2| it is shown that the connection between the physical energy distributions and the 



complex energy eigenvalues given by § p3 . 



2 Eigenvalue problems of the parabolic potential bar- 
rier 

2.1 The coordinate representation 

A simple and interesting model of an unstable system in quantum mechanics is the one- 
dimensional parabolic potential barrier. The Hamiltonian in quantum mechanics is 

where Vq G M is the maximum potential energy, ?7i > is the mass and 7 > is propor- 
tional to the square root of the curvature at x = 0. This Hamiltonian can be obtained 
from the Hamiltonian of the harmonic oscillator ( |A.1[ ) with the replacement 



u -^ Ti7- (2.2) 

Thus H is an essentially self-adjoint on a Schwartz space iS(M2:)- Also, two conditions are 
satisfied (see Appendix^. 



(i) 5(ffia) is an invariant subspace of H. 
(ii) H is continuous on 5(]Rk). 
The time-independent Schrodinger equation (|L^) is 

h^ CP , . (,. 1 2 2 



2m dx^ \ 2 

It is convenient to introduce the dimensionless variables 



u{x) + I Vo — -?Ti7 X u{x) = Eu{x). (2.3) 



Equation (|2.3|) now becomes 



i=l3x, I3=-J'-Y' (2.4) 

A.H(^. (2.5) 



^+(A + a« = 0. (2.6) 



Step 1. The asymptotic solutions 

We shall now obtain the asymptotic solutions of equation (p^). If we neglect A in ( |2.6| ) 
altogether, approximate solutions for large \^\ are 

u^{0 ~ ee^'^''\ (2.7) 



where n is an arbitrary number. Substituting (|2.7|) into (|2.6|) , we see that the values E^ 
of (^.51) associated with u^ are complex numbers: 



E+ e C_, E- e C+, 

where the symbols C± have the meaning 

C± = {z = x + iy\xeR, ye M±}, 
M+ = [0,+oo), M_ = (-00,0]. 

Since H is an essentially self-adjoint on ^(IRi:), m^ do not belong to a Lebesgue space 
L'^{M.x). For such eigenfunctions of H belonging to the complex energy eigenvalues, one 
can use the Gel'fand triplet which is formed by the Schwartz space iS(Ma;) as a nuclear 
space [|, §, 

5(M,)CL'(M.)C5(R,)\ (2.8) 



Let us treat the extension H^ of the Hamiltonian to the conjugate space iS(Mi.)^. We 
should be able to apply it to a generalized function u G ^(Mj:)^, the product H'^u being 
defined by [g ^ 

{v\H''u) = {Hv\u) 

for all functions v G 5(Ma;). Taking the a;-representatives, we get 



v{xy{H''u){x)dx = / {Hv){x)*u{x)dx. 

) J —OD 

We can transform the right-hand side by partial integration and get 

/■OO 

v{xy[H'^u){x)dx = / v{xy[Hu){x)dx, 



since w is a rapidly decreasing function, the contributions from the limits of integration 
vanish. This gives 

{v\H''u) = {v\Hu), 

showing that 

H^'u = Hu. 

Thus H^ operating to a generalized function has the meaning of H operating. 

The asymptotic solutions (|2.7|) will then be generalized functions or tempered distri- 
butions in iS(Ri.)^. We may verify that u^ G 5(]Ra.)^ by the Gauss- Fresnel integral, e.g. 

for any x^e-"'^'/^ ^ 5^^^) (m G Z+, a G M), 



00 



Xm^-oPx^/2\\n^±if5^xV2^^ = (m + n- 1)!!, 



27r 



(a^ ip i[] 



2\m+n+l 



when m -|- n is even and zero otherwise. 

Let us examine physical properties of the asymptotic states. The result of the mo- 
mentum operator p^ = —iHd/dx apphed to asymptotic solutions (|2.7|) is 

-in—u^{x) ^ ±}i(3^xu^{x). (2.9) 

dx 

Thus u~^{x) represents particles moving outward to the infinity x = oo, and u^{x) repre- 
sents particles moving inward to the origin x = 0. We shall see in § ^TT] that this space 
boundary condition is justified by the asymptotic behaviors of the probability currents. 

Step 2. The method of power series expansion 

We shall work out the whole solutions of equation ( |2.6| ) by the method of power series 
expansion from ( p.?!) . Put 

u±(e)=e±^«'/2iJ±(0, (2.10) 



introducing two new functions if^(^). Equation ( p.6| ) becomes 

^±2.«f^+(A*±.)«*-0. (2.11) 

We now look for a solution of these equations in the form of power seriesQ 

oo 

H^{0 = eY.''nC {4^0), (2.12) 

n=0 



in which values of the leading order s need not be integers. Substituting ( p.l2| ) in ( 2.11|) 
and picking out coefficients of C,"', we obtain 



sis - 1)4 = 0, 



s + l)scf = 0, } (2.13) 



{s + n + 2){s + n + l)c^+2 = M(2s + 2ri + 1) - A±] 






From the indicial equations of (2.13) 



s = or 1. (2.14) 

We can investigate the convergence of the series (p.l2|) on the same lines as the har- 
monic oscillator. The third of equations ( p. 13 ) give approximately, when n is large, 



cj n 



The series ( p. 12 ) will therefore converge like 

for large values of |^|. Thus, from ( p.lU| ), we have 

These results change the asymptotic behaviors in which u^{C,) goes off to infinity and m~(0 
arrives from infinity. There are therefore in general no permissible solutions of ( p^.3D for 
which the series ( p.l2| ) extend to infinity on the side of large n. The permissible solutions, 
which satisfy the space boundary condition, arise when the series (|2 . 1 2| ) terminate on the 
side of large n. The conditions for this termination of the series are 

cf =0 (2.15) 



^One cannot infer that u {^) have a definite parity in this step, even though the potential energy of 
(2.1) is a symmetrical of ^, since u'^iS.) are generally degenerate fi M. 



and 

A± = ^i(2s + 2m+l) (2.16) 

for some even m. In this case we get, from the recurrence formulas of (|2.13|), 



'^i ~ '^3 — ■ ■ ■ — c^_i — c^+i — c^+2 — • • • — U. 

Thus -ff^(0 cire the polynomials of degree s + m. Since the only possible values of s are 
or 1, we obtain from ( p.l6| ), by putting s + m = n, 

E^ = VoT^(n + ^^nJeC^ (n = 0, 1, 2, ■ ■ ■ ) ■ (2.17) 

This formulas give the complex energy eigenvalues for the parabolic potential barrier. 
For generalized eigenstates of H belonging to the generalized eigenvalues E^, the time- 
independent wave functions are 



.±^^, _ ^±^±.;3^.V2^±(^3,) e s{R,^\ (2.18) 



u„ix 



where H^{f3x) are the polynomials of degree n (see Appendix §), and B^ are the numer- 
ical coefficients. These numerical coefficients cannot be determined by the normalizing 
condition (|1.9|) , since a scalar product of generalized functions is not defined. In general, a 
generalized function is not necessarily normalized [0 . However we can determine B^ with 
the help of the following argument. Comparing above results with ([A.6|) and ( [A.7|) , we 



see that E^ and u^{x) are the same form as in the harmonic oscillator from the analytical 
continuation ( |2.2|) (see Appendix 0). From ( [A.8|) , we therefore determine that 



^" V(±2i)"n!^/S^J ^^'^^^ 

with the correction from our definition of H^{f3x) of Appendix p. 

Let us now see the properties of time-independent wave functions under a space in- 
version. Property (|B.3|) of Appendix shows that u^{x) are eigenstates of the parity. 



The validity of this result depends on the space boundary condition. 

The above results suggest that the eigenfunctions of the parabolic potential barrier 
are non- stationary states. It is obvious physically that the observational energy spectrum 
of this system is real and continuous. But the energy eigenvalues ( |2.17| ) are complex and 
discrete. The connection between the real continuous spectrum of observed energies and 
the complex discrete eigenvalues of mathematical eigenvalue problem will be dealt with 
in 5 0. 



2.2 The momentum representation 

The coordinate and the momentum representatives of a good function / G ^(Mx) are 
connected by the Fourier transformations [|l|] 






V2nh J-oo 

1 poo 



V27ih 






(2.20) 



Let us treat the extension T^ of the Fourier transformation to the conjugate space 
5(Ma;)^- We should be able to apply it to a generalized function u G ^(Mi;)^, the product 
J^^u being defined by [g ^ I 

(v\P''u) = {:f-'^v\u) 

for all functions v G 5(M^). Taking the representatives, we get 

/OD 
(^^^^v) {x)*u{x)dx. 

We can transform the right-hand side by the second of equations ( |2.20|) and get 

v{p)*(yP^u){p)dp = / v{p)*(yPu){p)dp, 



with the exchange of the order of the integrals by Fubini's theorem. This gives 

{v\P''u) = {v\Pu), 
showing that 

Thus JF^ operating to a generalized function has the meaning of JF operating. 



We note that the Fourier transformation and its extension applied to (|2.8|) give a new 
Gel'fand triplet 



5(M,.) C L2(M,) c S{ 



T 



T 



T'^ 



(2.21) 



cS(Rp) C L^ 



ip) c 5(Mp)^ 



which is a unitary equivalent to ( p^.Sj ) [|, 0, |] . Written in terms of p-representatives, the 
Hamiltonian ( |2.1| ) gives 



PhP~^ 



p" 
2m 



2fc2 



+ Vo + -mYh 



dp^ 



(2.22) 



an essentially self-adjoint on iS(Mp), and also the two conditions (|) and (0) of § ^TT] are 



satisfied. For generalized eigenstates of (|2.22| ) belonging to the generalized eigenvalues 
E^, the time- independent wave functions are 

<{p) =S„^e^^'^V/2^T(^p) e s{%)\ (2.23) 

5±=(-2)"| ^ ^ I , ^ = _^. (2.24) 

Note that the phase factor (— ^)" appears as a result of the Fourier transformations of 



3 Motion in the parabolic potential barrier 

3.1 The time evolution 

Our work in § |^ has been concerned with one instant of time. The present section will be 
devoted to the variation with time of the non-stationary states. 

We can get the time factors T^(t) corresponding to the complex energy eigenvalues 
E^ by substituting (|2.17| ) into (|1.6|) . They give 



where A^ are the numerical coefficients. The time-dependent wave functions ilj^{t,x) 
representing a non- stationary states of complex energy eigenvalues E^ are, according to 



PD 



V'±(t,x)=T±(t)M±(x) 

=A±B±e-*^«*/V('"+i/2)^*e±*'^'^''/2ff;^(/?x). (3.2) 



Before proceeding to discuss the physical meanings of the time factors (|3.1|) , we shall 
first verify the equation of continuity. From ( p.2| ) we see now that the probability densities 
are 

= l^n r \Bi\"e^^'-^'^^'H^iPx)H^iPx) (3.3) 

and the probability currents are 

j±(t,x) = ^ [i^i{t,xr {-zHdij^/dx) {t,x)]/m 

_ _L I /l±|2 I n±|2 zF(2n+l)7t 
—I— ■^^77 '^ n 

X 7 {xH^{(3x)H^{(3x) ± 2n(3-^^ [H^{(5x)HtAPx)\ ] . (3.4) 



Equations ( p.3| ) and ( p.4| ) depend on the time t which satisfy the equation of continuity 

^Pn(i,a;) + ^j±(t,x) = 0, (3.5) 

with the help of formulas (p.6| ) and ( p.7|) of Appendix p. Note that the numerical 
coefficients A^ and i?^ do not affect the validity of ( f3.5| ). 

We must examine how the probability currents jj^(t, x) behave for large values of |a;|. 
The second term in the {} brackets in ( |3.4| ) can be neglected as |x| -^ oo. We are left 
with 

These asymptotic behaviors of probability currents show that the index + (-0+, u^, etc.) 
means only outward moving particles and the index — {ip~, u~, etc.) means only inward 



moving particles. These asymptotic behaviors are the one assumed in § |2.1| for dealing 
with the space boundary condition. 

Let us now examine how the probability densities p^it, x) vary with time. From ( p.3|) , 
Pnit, x) will tend to infinity as t — > =Foo. The result is a necessary consequence of the form 
of time-dependent wave functions ( p.2| ). In any case the time-dependent wave functions 
ipnit,^) must not tend to infinity as t ^ =Foo, or they will represent states that have no 
physical meaning. The way of escape from the difficulty lies in the domains of the time 
factors ( p.l|) . We assume that T+(t) exists only on the future part, for which t > 0, and 
T~{t) exists only on the past part, for which t < 0, i.e. 

T+{t) =^+e-*^ot/?ie-(™+i/2)7t ^hen t >0, 
T-(t) =yi-e-*^ot/?ie(n+i/2)7t ^hen t <0. 

These equations can be combined into the single equations 

T±(t) = v4±e-*^«*/V("+^/2)^*^(±t), (3.6) 

where 9{t) is the Heaviside step function. We take the squares of the moduli of the time 
factors ( |3.6D 

|T±(t)|' = |A±|'e^(='"+^)^*^(±t). (3.7) 

To interpret the result (|3.7|) , we may suppose that the index + (■?/^+, T+, etc.) refer to the 
states of decay when t > 0, while the index — (■?/'", T~, etc.) refer to the states of growth 
when t < 0. This condition expresses mathematically that 

T+(t) G L^(M,), T~{t) e L2_(R,), (3.8) 



10 



where spaces L'^iM.t) are defined by 

LliRt) = {fe L\Rt) I fit) = a.e. t E R^}. (3.9) 

This condition will be referred to as the time boundary condition. 

These boundary conditions require that ip~{t,x) represents the growing state lying 
on the past part together with only inward moving particle and ip^it^x) represents the 
decaying state lying on the future part together with only outward moving particle. We 
must now obtain the probability of a transition taking place from state ■?/'" to state ip^ 
at time t = 0. This gives us the "selection rule" at time t = 0, and can be worked out 
from the following transition matrix elements 

Smn = {u:^\K)= U^{xyU~{x)dx, (3.10) 



the last definition expresses (m,^ | m„ ) in terms of ^-representatives. When we substitute 
for u^{x) their value given by (|2 18|) , equation ( p.lO|) gives 



^mn P 



/OD 
HiiiYH-XOe-'^'di 



with the help of ( p.4|) . The right-hand side becomes, from an application of formula ( p.8|) 
of Appendix 0, 

^mn ^ (^mn l3.iij 

with the help of ( p.l9| ). The transition matrix is then a unit matrix. This result expresses 
that the state ip~ (n = 0, 1, 2, ■ ■ ■ ) on the past part is connected with the state ip^ having 
the same label n on the future part. Thus our selection rule at time t = is that only 
those transitions can take place in which the complex energy changes from E~ to E^. 
Our problem now is to obtain the numerical coefficients A^ of (|3.6| ). We shall distin- 



guish between the growing state ipn ^^^ the decaying state Tp^ of the unstable system, 
with the growth and decay processes governed independently each by its own probability 
law []I[ . We shall require that the total number of growing particles is equal to the total 
number of decaying particles, i.e. 

p-{t,x)dt= / pl{t,x)dt. (3.12) 



By substituting for p^ here their values given by ( p.3|) , we now see that this condition 
( p.l2| ) is independent of x. Thus ( |3.12[ ) provides the condition for the time factor. From 
( P.12| ) we can require the normalizing condition for the time factor, as in equation ( |1.11|) 

of^ra 



/■oo 
|2 ,, / lm+/,x|2 



T-{t)\dt= \T^{t)\dt=l. (3.13) 



11 



If we now apply these conditions ( p.l3| ) to the time factors ( |3.6| ) , we have 



A^ = v/(2n+l)7, (3.14) 

where the phase factors are chosen unity.|^ 

Let us now apply the rules of ( |1.12|) for the interpretation of time factors. We find 
that the expectation values of t for the states ip^ are 



/oo 
t\T^it)\'dt 
-co 

/oo 
-oo 



-oo 

1 



(2n + l)7' 
respectively. Again, we find 



(3.15) 



/oo 
t'\T^itfdt 
-oo 

/oo 
-oo 



-oo 

2 



[(2n + l)7]2' 
Hence the positive square roots of the variance of t for the states ijj^ are 

= \/(*=>* - ({€f 

= 7 ^^- (3.16) 

(2n + 1)7 ^ ^ 

Thus (/^t)+ is equal to (^t)^, since the probability densities |T+(t)| and \T^{t)\ are 
symmetrical with respect to the center t = 0. We may call these quantities {^i}^ the 
mean lifetimes. Equations (|3.16|) inform us that the m,ean lifetimes are quantized in the 



unstable system with the Hamiltonian (2.1). These results are closely connected with 



Heisenberg's uncertainty relation between time and energy. This matter will be dealt 
with in the next section. 

Let us return to the time-dependent wave functions ( |3.2|) . We shall verify finally a 
time reversal. By application of property ( [B.4| ) of Appendix ^, we get 

^i{-t,xy = ^^{t,x), 



^We choose these phase factors so that the time-dependent wave functions are simple forms under a 
time reversal. 

12 



with the help of ( p.l4| ) and also of equation ( |2.19| ). Similarly we have in the p-representation 



since the phase factor in equation ( f^.24| ) and the parity of H^{j3p) cancel out. We see in 
this way that a time reversal occurs resulting in the interchange of states on the future 
part and on the past part. It is known as the reciprocity theorem H]. 



3.2 The energy distribution 

The time factor and the physical energy dependence of a function / G L^(Mt) are con- 
nected by the inverse Fourier transformations 



fiE)={:F-'f)iE) 
f{t)={Pf){t) 



V2^ 



V27rh 



f{E)e-'^''HE, 



(3.17) 



since the theory of relativity puts physical energy in the same relation to time as momen- 
tum to coordinate (see equations (p.20| )). 

We must now evaluate the physical energy dependence T^{E) for the time factors 
( p.6|) , and obtain the probability density of them. We make use of the Paley- Wiener 
theorem, which states that the inverse Fourier transformation is a unitary mapping from 



Lim 



onto H^, i.e. 



:F-'Ll{Rt) = H._ 



2 

±1 



(3.18) 



iJi denoting a Hardy space on the upper half-plane and H^ denoting a Hardy space on 



the lower half-plane. If we apply the first of equations ( |3.17| ) to the time factors ( p.6| ) and 
(Of), we then have 



f^E) 



At 



V2TThJ- 



'^''-^^^'^'/''e^^''+'/^^^'e{±t)dt, 



eH 



^/^E-Vo±i{n + 1/2)^7 " ^ ^' 
with the help of ( |3.14| ). Thus T^{E) has a simple pole at the point 

K)-«(r^ + 1/2)^7 eC_, 
and T~{E) has a simple pole at the symmetrical point 

Vo + i{n + l/2)H-feC+ 



(3.19) 



13 



with respect to the real axis. To obtain the original equations ( |3.6| ) we note that, from 
an application of the second of equations (p.l7|) to (|3.19|) , the integral being taken along 



a large semicircle in the lower half-plane when t > and in the upper half-plane when 
t< 0. 

By taking the squares of the moduli of ( |3.19| ) respectively, we obtain 

iT-ml'- iT+rnl'- ^ (n + l/2)/i7 

\TAE)\ -\TAE)\ --^E-Voy+[{n+l/2)njr ^^'^^^ 

These physical energy distributions are the Breit- Wigner resonance formulas, having the 
same resonance energy Vq and the half-widths 

r^ = (^ + ^) %• (3.21) 



Specially, the quantity F^ is what is sometimes called the decay width. Equations (|3.21|) 



inform us that the half-widths are quantized in the unstable system with the Hamiltonian 
( p.l|) . The integrals of ( p.20| ) with respect to E are 



2 



oo 



TAE)ydE= T+(E)rrfE = l, (3.22) 



corresponding to (|3.13|) , so T^{E) are normalized. 

Let us take the products of the mean lifetimes (Z\t)^ and the half-widths r^. From 
( P.16| ) and (|3.21| ) these products read 



m^r^ = I- (3.23) 

These results are independent of n. There are thus Heisenberg's uncertainty relations 
between the mean lifetimes and the half-widths. 

4 Discussion 

We have obtained some insight into the main features of the unstable system by making a 
study of the one- dimensional parabolic potential barrier. The results of solving the time- 
independent Schrodinger equation for the unstable system are complex energy eigenvalues 
and generalized eigenfunctions in Gel'fand triplet. These generalized eigenfunctions are 
non-stationary states. We have assumed that the time factor corresponding to a non- 
stationary state is a square-integrable function. This assumption requires that the states of 
complex energy eigenvalues in the lower half-plane exist only on the future part and those 
of complex energy eigenvalues in the upper half-plane exist only on the past part. They 
correspond to the states of decay and growth, respectively. The distributions of physical 
energy of these decaying and growing states are given by the Breit- Wigner resonance 
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formulas whose half-widths are determined by the complex energy eigenvalues. Thus 
the half-widths are quantized in the unstable system instead that the values of physical 
energy are quantized in the stable system (e.g. the harmonic oscillator, the hydrogen 
atom). These results are capable of comparison with experiments. The observational 
results will provide a confirmation of the complex energy eigenvalues and generalized 
eigenfunctions in Gel'fand triplet. 

We may discuss the comparison with the formulation given by Bohm and Gadella |0] . 
According to their formulation Gamow vectors (i.e. the states of resonance corresponding 
to a pair of simple poles of the S'-matrix) are represented by generalized functions of the 
following Gel'fand triplets, 



5±(Mi) C LliRt) C S, 



'±[i^t. 



T~ 



T- 



{T- 



A- 



C 



^1 



c 



z\+", 



where >S±(Mi) = 5(Mt) n L^(Mt) and Z\± = S{^e) n H\. The time evolution on A^"" 
is well defined only on the future part and that on Z\_^ is well defined only on the past 
part. We could, however, express the time factor as a Hilbert vector (§ ^, while we 
described the time- independent wave function by Gel'fand triplets (§ 0). We thus see 
that the origin of the separation of the domain of the time evolution lies in the property 
of the square-integrable of the time factor. 

In §§ |3.1| and p.2| we supposed that there is no interference between the growth and 
decay processes. In dealing with these processes, a further possibility to be taken into 
consideration is the direct scattering [|l|. In this case, we must first require the continuity 
condition for the time factor at time t = 

T-(0)=T+(0). (4.1) 

The reason for this condition ( |4.1|) is that the time-dependent Schrodinger equation (|1 . 1|) 
is linear in the operator d/dt. We can now require the normalizing condition for the time 
factor, instead of (|3.13|) 



/oo 
\T-it) + T^it)\'dt = l. 
-oo 



(4.2) 



This condition ([4.2|) means that the total probability of the state n being at any time is 
unity. The total energy dependence, namely T~{E) plus T^{E), is a real function of E. 
By taking the square of the modulus of the sum, we obtain 



\f-{E)+f+{E)\ 



[(n + 1/2)^7]= 



n{{E-Voy + [{n+l/2)n^]^y 



(4.3) 



Comparing this with ( |3.2CI| ), we see that the distribution ( [4. 3D has a narrow width (its 
quarter- width is equal to F^ given by ( |3.21|) ) and a sharp peak (its maximum value is 
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twice that of ( |3.20| )) at the same point E = Vq. This distribution is different from the 
Breit-Wigner resonance formula. It is necessary to re-examine experimental data of an 
absorption line induced by the scattering process. We will investigate how this difference 
would appear experimentally. 

The foregoing work shows how the properties of the time factor of the non- stationary 
state (of the unstable system) are similar to the properties of the time-independent wave 
function of the bound state (of the stable system), and vice versa. Further, the physical 
energy of the unstable system is the continuous spectrum running from — cxd to cxd, and 
the lowest energy does not exist. We should thus expect to find time operators for the 
unstable system in quantum mechanics. The subject for further investigation is to get 
these time operators and set up a relativistic quantum mechanics treating the time and 
coordinates as observables on the same footing. 

Appendices 

A The harmonic oscillator 

The Hamiltonian of the one-dimensional harmonic oscillator is 

H = -^^ + Vo + ^m^V, (A.l) 

2m dx'^ 2 

where Vq G M is the minimum potential energy, ?7i > is the mass and tu > is the 
angular frequency. This Hamiltonian is an essentially self-adjoint on a Schwartz space 
iS(Ma;). Further, two conditions are satisfied |0, 0. 

(i) 5(M^) is an invariant subspace of H. 

(ii) H is continuous on ^(IRj^). 

The time-independent Schrodinger equation (|TT^) is 

h^ <P ( \ \ 

7-^u{x) + Vo + -mJ^x^ u(x) = Euix). (A. 2) 

2m dx'' \ 2 J 

It is convenient to introduce the dimensionless variables 



Equation (|A.2|) now becomes 



mu . . ^, 

4 =0!X, a = W-^, (A.3) 

A.?(^-i«. (A.4) 



^+(A-a- = 0. (A.5) 
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Under the boundary condition u G 5(M^), this equation can be solved by the method of 
power series expansion. The energy eigenvalues are 

En = Vo+ (n + ^jnuem {n = 0,1,2, ■■■) (A.6) 

and the eigenf unctions are 

Unix) =iV„e-"'^'/2^„(«x) G S{R,), (A.7) 

1 

Nn = ( — ^^ ' . (A.8) 

Here Hn{ax) are the Hermite polynomials of degree n, and the constants A^„ are deter- 
mined by the normalizing condition ( |1.9| ). The set of eigenfunctions {tinj^o' Si^^n by 
( |A.7| ), forms a complete ortho normal system on a Lebesgue space Li^i^x) @|. 

The energy eigenvalues ( A.6 ) show that the probability densities of one-dimensional 



harmonic oscillator are time-independent, i.e. the eigenfunctions ( [A.7|) are stationary 
states. 

B Properties of the polynomials H^(^) 

The differential equations ( ^.11| ) with (2.16) 



^±2^^^T2^nH^ = 0, (B.l) 

where H^{^) are the polynomials of degree n. If we define them so that the highest power 
of ^ appears with the coefficient 2", the first few polynomials calculated from ( |2.13| ) are 

i7o±(0=l, i/i^(0 = 2e, Ht{0 = ^eT2t, 
HtiO =8^' T 12^e, ^4^(0 = 16^' T 48^e' - 12. 



These polynomials in general can be expressed as 



H^iO = iT^re^^^'^e^^^\ (B.2) 



Thus the parity of H^{C,) is 

HH-0 = i-TH^iO, (B.3) 

and the conjugate complex of H^{C,) is 

H^iO* = H^iO- (B.4) 
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An alternative way of defining the polynomials H^{C,) is as the generating functions 
S^i^^s) given by 



oo 



^±(e, s) = e^^(2^«-^') = Y^ ^^^iT^s^. (B.5) 



n=0 

From this definition we obtain the recurrence formulas 

H^^, - 2^H^ ± 2inHti = 0, (B.6) 

dH^/d^ = 2nHti. (B.7) 

Corresponding to the orthogonality relation of the Hermite polynomials, we now have, 
from (IB3D 



/oo 
HliiTH^iiy^'^'di = {±2iYn\V^m5mn, (B.8) 

-oo 

where 5mn is the Kronecker delta symbol. 
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